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ABSTRACT 



An analytical study was undertaken to determine two- 
dimensional wall conduction effects during film condensation 
on an internally finned, rotating heat pipe. An earlier 
Finite Element method was modified to solve the problem 
by using different numbers of axial increments and different 
numbers of elements. Resulting heat transfer rates were 
shown to be similar to those obtained earlier. 

A rotating heat pipe was tested using various copper 
condensers including two smooth-wall cylinders and a 0.5 
degree truncated cone. All condensers were tested at 
different rotational speeds using distilled water as the 
working fluid. In the case of the 1.46 inch diameter 
cylindrical condenser, data were also taken with Freon 113 
and ethanol. The heat transfer rate of each condenser was 
plotted against the saturation temperature of the vapor. 

Using the experimental results obtained in this thesis, 
together with the data of earlier theses, a comparison 
was made with the theoretical analysis of Roetzel. Agree- 
ment is reasonable, and Roetzel 's analysis can be used to 
predict overall rotating heat pipe performance to within 
± 20 percent for water. The data for Freon 113 and ethanol, 
however, fall approximately 35 percent higher than this 
theory. 
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I. 



INTRODUCTION 



A. THE ROTATING HEAT PIPE 

The rotating heat pipe is a closed rotating container 
capable of transferring large amounts of heat from rotating 
machinery. It consists of three main parts: (1) a cylin- 

drical evaporator, (2) a cylindrical or truncated-cone 
condenser and (3) a working fluid. The operation of a 
rotating heat pipe is shown schematically in Figure 1. 

When the rotating heat pipe is rotated above a critical 
speed, the working fluid forms an annulus in the evaporator. 
Heat input to the evaporator causes some of the working 
fluid to vaporize, and this vapor flows into the condenser 
transporting latent heat of vaporization. External cooling 
of the condenser allows the vapor to condense on the inner 
condenser wall. The centrifugal force due to the rotation 
has a component acting along the condenser wall which 
accelerates the condensate back to the evaporator. 

B. ROTATING HEAT PIPE LIMITATIONS 

Ballback [1] has studied the limitations imposed on 
the rotating heat pipe by various fluid dynamic mechanisms 
in a preliminary theoretical analysis of its operation. 

By using existing theories and experimental correlations, 
he was able to estimate the sonic limit, the boiling limit, 
the entrainment limit and the condensing limit of performance. 
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Figure 1. Schematic Drawing of a Rotating Heat Pipe 



A typical rotating heat pipe with the physical dimen- 
sions as shown in Table 1 was chosen in order to demon- 
strate the predicted performance. 

Table 1 

Specifications of a Typical Rotating Heat Pipe 



Length 


14.000 


inches 


Minimum diameter 


2.000 


inches 


Wall thickness 


0.125 


inches 


Internal half angle 


1.000 


degree 


Rotating speed 


2700 


RPM 



1 . Sonic Limit 

As the heat flux in the rotating heat pipe is 
increased, it is possible that the flow rate of the vapor 
down the container may be limited by the attainment of a 
choked flow condition for the vapor. This choking of the 
vapor flow will definitely limit the amount of energy the 
vapor can transport down the pipe, thus limiting the opera- 
tion of the pipe. The heat transfer mechanism employed in 
the rotating heat pipe depends on the latent energy of the 
working fluid. The heat transfer rate in the pipe may 
then be expressed as: 



m h - 
V fg 



( 1 ) 



where 



= mass flow rate of the vapor, Ibm/hr. 

From continuity, the mass flow rate of the vapor may be 
expressed by: 

m = p u A (2) 

V ^v V 

where 

u^ = velocity of the vapor, ft/sec 

2 

A = cross section area for the flow, ft . 



Then 

q. = p u A he 

^t '^v V fg 

The vapor velocity is taken as the sonic velocity, so 
that : 



u = c = /g kRT 
V ’o 



where 



(4) 



c = sonic velocity, ft/sec 
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g = acceleration of gravity, 32.1739 ^=- 7-5 9 

^o '3 j > Ibf-sec'^ 

k = the ratio of specific heats, 

R = gas constant, ft-lbf/lbm-deg R 

T = absolute temperature, deg R 

Using equations (3) and (4), the sonic limit of a 



typical rotation heat pipe can be computed, as shown in 
Figure 2. 

2 . Boiling Limit 

Kutateladze [2] has postulated that the transition 
from nucleate boiling to film boiling is purely a hydro- 
dynamic phenomenon. Based on this assumption he has 
obtained the correlation formula for predicting the burnout 
heat flux: 



1/4 



q. = K/^ a. h- 

^t V b fg 



{ag[p^ - P^] } 



(5) 



where 



A, 



b 



2 

heat transfer area in the boiler, ft 



P 



density of the vapor, Ibm/ft 



3 



V 




a 



surface tension, Ibf/ft 
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Heat Transfer Rate, q BTU/hr 



1.7 



17.2 



2 

Saturation Pressure, Ibf/in 

2.9 4.7 7.5 11.5 




Figure 2. Operating Limits of a Typical, Water-Filled, 
Rotating Heat Pipe 
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r a.( 



g 



p 



£ 

j- 



2 

acceleration of gravity, ft/hr 

3 

density of fluid, Ibm/ft 



K 



constant value. 



Kutateladze ' s experimental data indicate a value of 



K between 0.13 and 0.19. Using equation (5), the boiling 
limit of a typical rotating heat pipe can be computed, as 
shown in Figure 2. 

3 . Entrainment Limit 

Sakhuja [3] studied the flooding constraint in 
wickless heat pipes, where he assumed that: 

a) the vapor and liquid are in counterflow motion 
and are exposed to each other; 

b) the falling film of liquid (under gravity) is 
subjected to a shear stress by counter flowing 
vapor, which tends to retard the film; 

c) the liquid film remains smooth and stable, and 
the shear stress is usually small; 

d) the buoyancy forces, due to density differences 
between the vapor and the liquid, are responsible 
for maintaining the counter flow; and 

e) the buoyancy forces are balanced with dissipative 
effects which are proportional to momentum 
fluxes of vapor and liquid streams. 

He has obtained a correlation for flooding: 




2 

A C h^ /gD ( p,r - p ) p 

X f g ^ '^f '^v 

[ 1 + 



( 6 ) 
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where 



= heat transfer rate, BTU/hr 

2 

A = flow area, ft 

X 

C = dimensionless constant, 0.725 for tube 
with sharp edged flange 

D = inside diameter of the heat pipe, ft 

2 

g = acceleration due to gravity, ft/hr 
= density of the fluid, Ibm/ft^ 

3 

= density of the vapor, Ibm/ft 

hgg = latent heat of vaporization, BTU/lbm. 

Using equation (6) , the entrainment limit of a 
typical rotating heat pipe can be computed and is shown 
in Figure 2. 

4 . Condensing Limit 

The solution for condensation in the rotating heat 
pipe has been carried out by Ballback [1]. He modeled the 
condenser section of the heat pipe as a rotating, truncated 
cone, and an equation to compute the condensing limit was 
derived: 









2 p? h [T -T 
2 f f fq s w j 



Ug sin (|) 



{ [R + L sind)] - R 
o c 



8/3 8/3 3/4 

} (7) 



where 

total heat transfer rate, BTU/hr 

thermal conductivity of the condensate 
film, BTU/hr-f t-deg F 
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3 

= density of fluid, Ibm/ft 

0 ) = angular velocity, 1/hr 

= latent heat of varporization , BTU/lbm 

Tg = saturation temperature, deg F 

= inside wall temperature, deg F 

= viscosity of fluid, Ibm/ft-hr 

= half cone angle, degrees 

R = minimum wall radius in the condenser 
° section, ft 

L = length along the wall of the condenser, ft. 



Equation (7) is a function of the physical proper- 
ties of the working fluid, the geometry of the rotating 
heat pipe and the operating speed. Using this equation, 
the condensing limit for a typical rotating heat pipe can 
be computed. This result is also shown in Figure 2. 

It is clear from the results of Figure 2, that for 
the geometry and working fluid chosen, and under normal 
circumstances, the heat transfer rate is limited primarily 
by the amount of heat which can be transferred through the 
condenser wall. However, other limitations may become 
important as heat pipe geometry and operating conditions 
are altered. 



C. THEORETICAL INVESTIGATIONS OF FILM CONDENSATION 
ON SMOOTH ROTATING CONDENSER SURFACES 

Sparrow and Gregg [4] were the first to study condensa- 
tion on a rotating surface. They studied condensation on a 
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rotating disk at uniform surface temperature, T^, situated 
in a large quiescent body of pure saturated vapor. They 
assumed that the condensed liquid formed a continuous film 
on the disk and that fluid in this film moved radially 
outward along the disk under the action of the centrifugal 
force field. They arrived at the following solution for 
the situation where energy convection and acceleration terms 
are negligible: 



w 



= 0.904 {- 



Pr 
AT/h 



fg 



(8) 



where 



h 

k 



f 



V 

to 

Pr 

y 

C 

P 

^w 

"^sat 

AT 



fg 



local heat transfer coefficient, 
BTU/hr-ft^-deg F 

thermal conductivity of condensate, 

BTU/hr-f t-deg F 

2 

kinematic viscosity of condensate, ft /hr 
angular velocity, 1/hr 
Prandtl number, yC /k^ 

absolute viscosity of condensate, Ibm/ft-hr 
specific heat of condensate, BTU/lbm-deg F 
disk surface temperature, deg F 
saturation temperature of vapor, deg F 
T . “ T , deg F 

latent heat of vaporization, BTU/lbm. 
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Equation (8) predicts that the local heat transfer 
coefficient is constant across the disk and does not vary 
with radial position. 

Sparrow and Hartnett [5] continued the research 
using a rotating cone with the same assumptions used earlier. 
Their result for the heat transfer coefficient was: 



h 



cone^^disk 



1/2 

(sin <p) 



(9) 



where ()> is the internal half -angle of the cone. As pointed 
out by the authors, however, equation (9) applies only to 
cones that are not too slender. Using their result, the 
mean Nusselt number for a cone can be written as: 



Nu 



m 



hL 

k 



= = 0.904 ( 



2.2 4 
w sin (p L 

2 

V 



Pr 



1/4 



C AT/h. ' 
P fg 



( 10 ) 



Leppert and Nimmo [6,7] studied laminar film 
condensation on a surface normal to body or inertial forces 
such as a finite horizontal surface, or a rotating drum. 

They found that the local heat transfer coefficient varies 
with position, and the mean Nusselt number, for a constant 
surface temperature, with negligible surface tension effects, 
is : 



Nu 

m 



1/5 



0.64 Sh 



( 11 ) 



ijiloa 






onciai 



where 



Nu = mean Nusselt number. 



h L 
m 



m 



Sh = Sherwood number. 



5 

k AT y 



Equation (11) can be applied to the case of 

condensation on the inside of a rotating cylinder of radius, 

2 

r, where g is replaced by rw . 

The first theoretical investigation conducted at 
the Naval Postgraduate School was performed by Ballback 
[1] . His analysis was a Nusselt type analysis for film 
condensation on the inside surface of a rotating truncated 
cone of half angle cone, 4>, with a geometry defined in 
Figure 3. His assumptions were similar to those used in 
the Nusselt analysis of film condensation on a vertical 
wall, namely: 

1) steady state condition, 

2) film condensation, not dropwise condensation, 

3) laminar flow of condensate, 

4) constant fluid properties, 

5) no condensate subcooling, 

6) no momentum change through the condensate, 

7) no liquid-vapor inter facial shear stresses, 

8) linear distribution of temperature through the 
condensate film, 

9) negligible force of gravity, and 

10) constant wall temperature. 
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Additionally Ballback assumed: 



11) 


a very small film thickness, 6 (x) , 




12) 


no radial distribution of pressure 
vapor space. 


in the 


13) 


negligible density of the vapor. 




14) 


negligible velocity gradients in the 
circumferential direction, and 


15) 


very small slope of the condensate 
d5/dx when compared with tan (j). 


film. 



He used the boundary condition, that 5 (x) must be 
zero at X = 0, and arrived at the following solution for 
the local film thickness: 



6(x) = 



3 ’"f ^f ^"^s ""^w^ 

{3 f f . s ^[ 1 - { 



R 



■7 “2 — 2 — ^ 

p - w sin <p h 



fg 



R + X sin d) 
o ^ 



8/3 1/4 

} ]} 



( 12 ) 



where 



6 (x) 

k. 



0) 

T 



w 



fg 



film thickness, ft 

thermal conductivity of the fluid, 
BTU/hr-ft-deg F 

viscosity of the fluid, Ibm/ft-hr 

density of fluid, Ibm/ft^ 

angular velocity, 1/hr 

inside wall temperature, deg F 

saturation temperature, deg F 

latent heat of vaporization, BTU/lbm 

half cone angle, degrees 

minimum wall radius in the condenser 
section, ft. 
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The local heat transfer coefficient, h(x), is 



given by: 



^f 


2 sin 4. 1/4 

3 ,, , "o 

+ X sin 4)^ ^ 

o ^ 

(13) 



The mean heat transfer coefficient over the whole 
surface area of the truncated cone is given by: 



h 

m 


1 ^ 

i.L(2R^ + Lsln 4) / 2iih(x) (R^ + x sin 4.) dx 

(14) 


where 





L = length along the wall of the condenser, ft. 

Integrating equation (14) yields a value for the 
mean Nusselt number: 



h L 

Nu - 

m k 


0)2 l2 r2 1/4 

(15) 



where 



G(S) 


= {(1 + 6)®/^-!}^/“ ^ 
/F(2 + B) 
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and 



6 



L sin 
R 



o 



4 > 



Equation (15) can be put into the form: 



Nu 

m 






G(g) 

/B 



(17) 



where (Nu ) „„ is the result of Sparrow and Hartnett [5] 
m on 

for a cone, as given by equation (10) . The above result 
was confirmed as a special case of the later work of Dhir 
and Lienhard [8] , who studied laminar film condensation on 
axisymmetric bodies in nonuniform gravity fields. 

Daley [9] was able to improve on the theory intro- 
duced by Ballback by consideration of the thermal resis- 
tances in the condenser wall and on the outside surface. 
Both Ballback and Daley, however, neglected the drag due 
to vapor friction on the counter flowing condensate film. 

Newton [10] extended the analysis performed by 
Ballback and Daley to include the effects of vapor pressure 
drop and interfacial shear stress between vapor and conden- 
sate. Using the same geometry in Figure 3, he assumed that 
the momentum changes within the condensate film were small 
and arrived at the following differential equations for 
condensate film thickness, vapor velocity and heat transfer 
rate : 
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dx 



2 (X) P»(x) M.R(x) 



tan <p - 



2 2p v(x) 

P^ (X) COS4) [pgU) R(x) + — ] 



(18) 



dv 

dx 



2 (R + xsin 4) (T - T ) 
o s 



„ f6 (x) , t . Ill, 

P L— 1 + 1 — + ttJ n - 

h fg 



2 

R(x) 



2v (x) 
R(x) 



[sin 4) 



cos 



4 > 



dx 



] 



(19) 



^ = 
dx 



2it (R + xsin 4> ) (T - T ) 
o ^ s “ 

~<S (x) ~ t 7”i 

k h 
f w 



( 20 ) 



where 



Pj^(x) 

P^M 



V (x) 
T 



t 

6 (x) 
k^ 



<S^x) 2 4 

(R + xsin tp )— ^ — - -^cos 4> <S (x) , 

^ 3 

(R + xsin 4> ) - - cos p 6 (x) , 

o 2 — 2 

velocity of vapor, ft/hr 

saturation temperature of vapor, deg F 

ambient coolant temperature, deg F 

thickness of the condenser wall, ft 

film thickness of the condensate, ft 

thermal conductivity of condensate film, 
BTU/hr-ft-deg F 
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R(x) 

<P 



V 



thermal conductivity of condenser wall, 
BTU/hr-ft-deg F 

external heat transfer coefficient, 
BTU/hr-ft^-deg F 

latent heat of vaporization, BTU/lbm 
friction factor, dimensionless, 
vapor Reynolds number, 2R (x) p^v (x) 

radius along the condenser wall, 

R + xsin (p - 5cos (J) 
o 

half cone angle, degrees 

3 

density of vapor, Ibm/ft 
viscosity of vapor, Ibm/ft-hr. 



The integrations were started by assuming initial 

conditions 6 = 6 . and v = 0 at x = 0. With these initial 

1 

d 6 

conditions - 5 — = tan d) at x = 0. Newton obtained results 
dx ^ 

for half cone angles of 0, 0.1, 0.2 and 0.3 degrees. As 
the half cone angle increased, the solution of the differ- 
ential equations became extremely sensitive to the initial 

values of 6 . . He also noted that his results for 6 (x) 

1 

showed that (except near the ends of the condenser) 
d 5 

^ << tan p. He, therefore, made this approximation and 
simplified the model further using an order of magnitude 
analysis of the terms in equations (18) and (19) . With 
these simplifications, equations (18) and (19) were 
reduced to: 



30 







■^v T: 











;■*: 


£ . i£tt 



m 



- - 


m 

■^od| 




hi 


rj 


' la 


«* 




i.- 


' .'iRfll 







i: f 






I 



dv 

dx 



2 (R + xsin (j)) (T - T ) 
o s °° 

p + ir~ rIx) 

k h fg 
f w ^ 



2v(x)sin (p 
R(x) 



( 21 ) 



3 

<S (x) 



2 R(x)sin R(x) + 



2p v(x) 2 , 

-^1 - p^v(x)flR(x)ii|l-l 



y^R(x) 



= 0 



( 22 ) 



Equation (21) was numerically integrated using the 
initial condition that 5^^ = 0 and v = 0 at x = 0 and using 
the cubic equation for 6 (x) , equation (22). After the 
integration was completed for 5 (x) and v(x), the energy 
equation, equation (20) , was integrated over the length of 
the condenser to find the total heat transfer rate of the 
heat pipe. This theory was compared to the experimental 
results of Schafer [15] and Tucker [20] , predicting perfor- 
mance to within ±20 percent for water, ethyl alcohol and 
Freon 113. 

Daniels and Al-Jiamaily [11,12] extended the analysis 
of Ballback by including drag effects due to the vapor 
friction. Their solution for the local film thickness 
distribution along the condenser is given by: 



sh . hr - he . 

X X 3 X X 4vx x 



1=0 



(23) 
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where 



Sh 

X 



Dr 

X 



Re 



vx 



P 

6 (x) 



0 ) 



u 



g 

h 



fg 



Sherwood number, 



p (oir-g) sincjih^ x 

£2 

4ykAT 



px h_ X cos 

drag friction number, — ? ■ 



ykAT 



3 






pu X cos (j) 

two phase Reynolds number, — ~ — 

liquid density, Ibm/ft^ 

film thickness, ft 

angular velocity, 1/hr 

half cone angle, degrees 

shear stress at liquid vapor interface, 
Ibf/ft^ 



liquid condensate velocity, ft/hr 

2 

acceleration of gravity, ft/hr 
latent heat of vaporization, BTU/lbm. 



At any section x in the condenser, since heat is 
transferred through the condensate layer by conduction 
only, the local heat flux = kAT/6 (x) and equation (23) 
can be re-arranged as follows : 



+ jRe (M1)2q2 ^ 1q^ _ Sh = 0 (24) 

X 4 vx X X 3 X X X XX 
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where 



Q = local heat flux, BTU/hr-ft^. 



This equation relates the three main parameters 

controlling the rotating heat pipe operation assuming 

fixed geometry and working fluid. The mean Nusselt number, 

Nu , for this analysis was, 
m 



Nu 

m 



4 

3 




L 

L 







(25) 



where suffix L refers to quantities evaluated at the end 
of the condenser section x = L. 

Marto [13] compared the heat transfer rates of 
Daniels and Al-Jumaily with Ballback and the approximate 
numerical solution to equations (20) , (21) and (22) . For 

the case of Freon 113, equations (20) , (21) and (22) , 

which include both effects of interfacial shear as well as 
vapor pressure drop, gave slightly lower heat transfer 
rates than the analysis of Daniels and Al-Jumaily which 
includes only interfacial effects. Also the analysis of 
Ballback predicts slightly larger results than the computer 
analysis but lower results than those of Daniels and 
Al-Jumaily. This may be due to differences in the Freon 113 
fluid properties used or due to the geometrical inconsistency 
in matching all the geometrical dimensions as stated in the 
analysis of Daniels and Al-Jumaily. For the case of water. 
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Ballback's results again predict a slightly higher heat 
transfer rate than the approximate numerical solution, due 
to his neglecting effects of interfacial shear and vapor 
pressure drop. When tjiese effects are neglected in the 
approximate numerical solution, the results agree very 
well with those of Ballback. For these reasons, the 
effects of interfacial shear and vapor pressure drop can 
safely be neglected. 

Roetzel [14] studied condensation in steam heated 
fast rotating paper drying drums, as shown in Figure 4. 

In his analysis, the heat transfer is determined mainly by 
the method of condensate removal. The condensate is collected 
in a circumferential groove on the inside surface, and from 
there it is removed by means of a stationary syphon with 
its tip in the groove. The mean film thickness is slightly 
smaller than with a rotating syphon because the condensate 
is sucked off from a lower level. The solution for the 
case of a drum with an infinite radius of curvature (The 
inside surface is exactly cylindrical.) with a uniform 
temperature difference is given by: 

1/5 

Nu = (0.86 E) (26) 

00 



where 



Nu^ = Nusselt number for the cylindrical drum. 
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b ph L 

"h = 6^ Jj - 

f m 

thermal conductivity of condensate, 
BTU/hr-ft-deg F 



n 



Sh 

b 

P 

h 



fg 



V 



= Sherwood number 

2 

= acceleration normal to the wall, ft/hr 

3 

= density of condensate, Ibm/ft 
= latent heat of vaporization, BTU/lbm 
= total length of condensate, ft 

2 

= kinematic viscosity of condensate, ft /hr. 



This equation can be written in terms of the 
Sherwood number : 



Nu = 0.678 Sh 

oo 



1/5 



(27) 



and compares favorably to the analysis of Nimmo and Leppert, 
equation (11) . 

For the case of finite radius of curvature, Roetzel's 
solution becomes: 



„ ,,2Le/3/16 



(0 



36 ^) 23 / 20 , 4/23 



(28) 
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where 



R = radius of curvature, ft, and 
L = total flow length of condenser, ft. 

This equation can be put in terms of Sh, Nu^o, R, and L as: 

Nu = Nu [1 + 1. 924 . (29) 

Thus the Nusselt number can be directly determined when 
the ratio L/R and the value of the dimensionless group, E, 
are known. For the truncated cone, rotating heat pipe, as 
shown in Figure 3, equation (29) can be applied to compute 
the heat transfer by defining the radius of curvature 
R = L/tan <() ; the solution then becomes: 

Nu = Nu [1 + 1.924 tan <f) j^23/80 ]4/23 

m oo '• 



where 



(}) = half cone angle, degrees 

Nu^ = Nusselt number of the cylindrical drum, 
equation (11) . 

Equation (30) has never been compared with experi- 
mental data, so using the existing rotating heat pipe data 
together with new experimental data at different half cone 
angles will provide a check on its usefulness. 
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These theoretical investigations are summarized 
below: 

1. The vapor pressure drop and inter facial shear 
stress effects between vapor and condensate can 
be neglected. 

2. The analysis of Roetzel [14] appears to be the 
most appropriate for use in a rotating heat pipe. 

3. The important parameters are: 

(a) speed, 

(b) geometry, and 

(c) Sherwood number, which includes the properties 
of the working fluid. 

D. ANALYSIS OF FINNED ROTATING CONDENSER SURFACES 

The use of axial fins on the inner surface of the con- 
denser can improve heat transfer performance through two 
mechanisms; a decrease in condensate film thickness near 
the apex of each fin and an increase in inner wall surface 
area. Schafer [15] developed an analytical model of the 
internally finned condenser for the case of a fin of tri- 
angular profile as shown in Figure 5, with a specific 
geometry as shown in Figure 6. In the analysis, he assumed 
one-dimensional heat conduction from the inside to the 
outside surface and a linearly decreasing surface tempera- 
ture profile from the fin apex to the midpoint of the 
trough. Results from his analysis indicated significantly 
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Figure 5, Internally Finned Condenser Geometry 

Showing Fins, Troughs and Lines of Symmetry 
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Figure 6. Condenser Geometry Considered in Analysis 
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improved heat transfer rates for the finned condenser 
compared to the smooth condenser. 

Corley [16] developed a numerical solution for this 
case using the Finite Element method to account for two- 
dimensional conduction effects within the wall for the 
triangular fin profile and a parabolic wall surface tempera- 
ture distribution along the fins. His results indicated 
significantly improved heat transfer performance with heat 
transfer rates about 75 percent above those predicted by 
the one-dimensional linear temperature distribution analy- 
sis of Schafer [15] . However, Corley noted that an error 
of 50 percent in the calculation of apex heat transfer rate 
might then give a total heat transfer error of as high as 
15 percent. Two possible solutions to the boundary con- 
dition problem were recommended. The first was to increase 
the number of Finite Elements in the analysis, and the 
second solution was to consider a wall geometry as shown 
in Figure 7 . 

E. THESIS OBJECTIVES 

The objectives of this thesis were therefore: 

1) to find out which existing theory best correlates 
Rotating Heat Pipe data by comparing these theories 
with existing data and new experimental data of 
this thesis, and 

2) to further study finned condensers by improving 
the Finite Element results obtained earlier. 
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Figure 7. Alternate Condenser Wall Geometry 
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II. FINITE ELEMENT ANALYSIS AND RESULTS 



A. REVIEW OF EARLIER ANALYSES 

Schafer [15] studied the one-dimensional heat transfer 
solution, and Corley [16] studied the two-dimensional heat 
transfer solution to finned condensers. Both used the 
same assumptions and boundary conditions, based upon the 
work of Ballback [1] , which are similar to those used in 
the Nusselt analysis of film condensation on a vertical 
wall. The more important of these assumptions are listed 
below: 

1) steady state operation, 

2) film condensation, as opposed to dropwise condensation, 

3) laminar condensate flow along both the fin and 
the trough, 

4) static balance of forces within the condensate, 

5) one-dimensional conduction heat transfer through 
the film thickness (no convective heat transer) , 

6) no liquid-vapor interfacial shear stresses, 

7) no condensate subcooling, 

8) zero heat flux boundary conditions on both sides 
of the wall section, as shown in Figure 6, 

9) saturation temperature at the fin apex, 

10) zero film thickness at the fin apex, 

11) negligible curvature of the wall. 

The complete development of Schafer's finned condenser 
heat transfer theory is presented in his thesis [15] . His 
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results are used for comparison purposes in this 
thesis . 

Corley's study [16] considered an internal fin geometry 
as shown in Figure 6. Two-dimensional conduction through 
the condenser wall was considered in the heat transfer 
analysis by using a Finite Element solution. Instead of 
using a linearly decreasing temperature profile along the 
wall from the apex of the fin to the midpoint of the 
channel as Schafer did, a parabolic distribution was assumed 
to calculate average nodal convective coefficients. The 
eight node finite element mesh used in the two-dimensional 
analysis of the axially finned condenser is shown in Figure 
8. A computer program for the case of two-dimensional con- 
duction was developed using a modified version of the three 
dimensional, isoparametric Finite Element program assembled 
by Lew [17] . 

A parabolic temperature distribution is assumed to exist 
along the fin surface using the boundary condition, T = T 

s 

at z = 0, since 5(0) = 0. Therefore; 



T (z) 
w 




s 



or 




2 




( 31 ) 
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Figure 8. Condenser Geometry Considered in 

Two-Dimensional Analysis with 2 Finite 
Elements 
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where 



and are unknown constants, and 

z is the distance along fin surface from apex of 
fin, ft. 

From Schafer's analysis [15] a differential equation 
for fin condensate thickness, 6(z), can be derived; 



6^(z)d6(z) = 


k^ (T^ - T (z) ) u, dz 
f s w £ ,32, 

2 2 

(R + xsin (t)h^ cos d) cos a 
f o fg 



where 



6 (z) 


the fin condensate thickness, ft 


11 

M-l 


thermal conductivity of working fluid, 
BTU/hr-ft-deg F 


Pf 


3 

density of working fluid, Ibm/ft 


11 

4-1 


viscosity of working fluid, Ibm/ft-hr 


0) = 


angular velocity, 1/hr 


11 

0 


minimum radius, ft 


X = 


distance along apex of fin, ft 


Z = 


distance along the fin, from the apex, ft 


(j) 


half cone angle, degrees 


a = 


fin half -angle, degrees 


fg 


latent heat of vaporization, BTU/lbm. 
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Equation (31) is then substituted in (32) , and the 
result is integrated from z = 0 to z to give an equation 
for 6 (z ) ; 



<S(z) = ( 



4Kf(-a^ Vf 



2 2 

P-O) (R + xsin 4>)h- cos 4» cos a 
to rg 



-] 



1/4 



(33) 



The local convective heat transfer coefficient assuming 
one-dimensional heat conduction through the thin condensate 
film is then: 



h(z) 



f in 



0 (z) 



3 2 2 

k-p^o) (R +xsin4>)h- cos <t> cos a 

2 _£2 , ] 



1/4 



4pj(-ai - 3 - - 



1 2 



-) 



(34) 



which is applicable from z=0toz*=z^- 6*/cos a , 
where 6* is the trough condensate thickness, as shown in 
Figure 8. Across the width of trough, from z = z^ to the 
mid channel position, the heat transfer coefficient is: 



trough 6 * (x) 



(35) 



where 6*(x) is considered constant across the trough width. 
Based upon a check of the wall temperatures generated from 
Schafer's computer program, the fluid temperature difference 
across the width of the trough was predicted to be small. 
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Fluid thermal conductivity values were assumed to remain 



relatively constant. As a result, is assumed 

to be constant across the trough width and is calculated 
using equation (35) and a representative average inner 
wall temperature to determine k^. The convective coeffi- 
cient along the outside surface of the wall, h is also 

assumed to be constant. The modified trough mass flow 
equation becomes : 






2 2 

p 0 ) (R^ + xsin 4>) 5 * (x) sin (}> 

TTT 



•[6*(x)e +5*(x)tan a] 

(36) 



Within the computer program, the heat transfer rate 
for the entire condenser is determined by dividing its 
axial length into one hundred equal length increments of 
length. Ax (See Figure 6.), solving for the incremental 
heat transfer rate using an iterative procedure, and 
summing these incremental rates to yield a total heat 
transfer rate. A flow chart description of the computer 
program used in Corley's analysis is given in Figure 9. 

At each increment, constant values of a^^ and b^ are deter- 
mined using a Lagrangian interpolation fit of the tempera- 
tures, at nodes 1, 3 and 6 as calculated from the previous 
iteration. For the first increment, a^ and b^ are calcu- 
lated from an initial guess of the temperatures at these 
nodes. Convective coefficients at several locations along 
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Figure 9. Computer Program Flowchart of Two-Dimensional 
Conduction Analysis of Corley [16] 
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the fin are determined using equation (34) . Average 



nodal values, and the average h of the 

left-hand boundary region associated with node 6, in Figure 
8, are calculated using the position dependent values of 
h(z)^^^ in a Simpson's Rule approximation. In the trough, 
is calculated from equation (35) using the 5* 
value calculated from the previous increment. This value 
is used as the average convective coefficient for nodes 9 
and 12 and the right hand boundary region of node 6. The 
average convective coefficients calculated for the left- 
and right-hand boundary regions of node 6 are weighted 
by their respective region lengths. The two values are 
summed to yield the average convective coefficient for node 
6, The remaining boundary conditions are applied, 

and the incremental heat transfer rate per unit condenser 
length, , and nodal temperature values are generated 
from the finite element subroutines. If the resultant 
heat transfer rate converges to within the criterion of 
0.05 percent change from the value of the previous iteration. 




<< 0.0005 , 



(37) 



then the incremental problem is deemed solved and 0^ is 
set equal to the final iteration, . If, on the other 
hand, the convergence criteria is not met, newly calculated 
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nodal temperature values are used to establish a new set 
of inner wall convective coefficients, and the finite element 
solution is repeated until the heat transfer rate converges 
to within the acceptable 0.05 percent tolerance. Because 
Qj is initialized at zero for each increment, the iteration 
procedure is repeated at least twice for every increment. 
Incremental mass flow rate is calculated using: 



AM 



tot 




(38) 



This value is added to the mass flow rate from previous 
increments. The subsequent interval's trough condensate 
thickness, 6*(x), is calculated using a polynomial root- 
finder subroutine to solve equation (36) . This incremen- 
tation is continued until the entire length of the conden- 
ser is traversed. Incremental heat transfer rates are 
summed, and a total heat transfer rate becomes: 



'tot 



= 2 N 



fin 



100 

z 

i=l 



(Q^- Ax) 



(39) 



where number of fins along the inner wall 

circumference. To start the incrementation, two initial 
assumptions are made. The initial value of 6* = 6| is 
taken from the analysis by Sparrow and Gregg of condensa- 
tion on a rotating disk [4] . In addition, at the first 
increment, initial guesses of fin temperatures at nodes 1, 
3 and 6 are made in order to begin the solution. 
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B. MODIFICATION TO THE FINITE ELEMENT ANALYSIS 



This thesis is a continuation of Corley's [16] study 
of the addition of internal axial fins to the condenser 
surface of a rotating heat pipe. He recommended that one 
way to approximate more closely the actual heat transfer 
rates would be by increasing the number of elements for 
the same geometry. Since Corley used one hundred equal 
increments of length to compute the heat transfer rates, 
the first step in this thesis was to reduce the number of 
increments to fifty and twenty five, to minimize the com- 
puter time. Using Corley's program, at a speed of 3000 

RPM, with external convective coefficients of 1000, 1500 

2 

and 5000 BTU/hr-ft -deg F and fin half angles of 10, 30 
and 50 degrees, the results for water as the working 
fluid are shown in Tables 2, 3 and 4. It is clear that 
reducing the number of increments to twenty five had little 
effect upon the heat transfer results (approximately 3 
percent error) . 

To reduce any error in the calculation of the heat 
transfer rate at the apex of the fin, the number of Finite 
Elements has to be increased between the apex of the fin 
and the trough, as shown in Figure 10. The geometry becomes 
3 elements and 18 nodes. The computer program was therefore 
modified to solve this case, using the same boundary condi- 
tions and assumptions. These results are shown in Tables 5 
through 10. Upon dividing element in Figure 10 into 
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TABLE 2 



Theoretical Finned Condenser Heat Transfer 
Rate Versus Number of Axial Increments 



h ^ = 1000 
out 


BTU/hr-ft^-deg 


F 


At 3000 


RPM 


Fin Half 


No. of T . 


No. ( 


of Axial 


Increments 


Angle 


• 5 

Fins 


25 


50 


100 


10 


276 100 


10273.1 


10039.0 


9922.5 




150 


27355.8 


26732.6 


26422.4 




200 


44427.2 


43415.0 


42911.2 


30 


84 100 


10095.3 


9865.3 


9750.8 




150 


26866.0 


26253.8 


25949.1 




200 


43622.0 


42627. 9 


42133.2 


50 


40 100 


9870.9 


9643.5 


9530.3 




150 


26255.7 


25650.8 


25349.9 




200 


42623.3 


41641.2 


41152.6 
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TABLE 3 



Theoretical Finned Condenser Heat Transfer 
Rate Versus Number of Axial Increments 

= 1500 BTU/hr-ft^-deg F At 3000 RPM 



Fin Half No. of T , No. of Axial Increments 

^ 1 Ti- sat cr, inn 



Angle 


Fins 


sac 


25 


50 


100 


10 


276 


100 


15195.2 


14849.2 


14677.1 






150 


40427.4 


39507.2 


39049.2 






200 


65632.8 


64138.9 


63395.3 


30 


84 


100 


14808.4 


14471.6 


14304.0 






150 


39368.6 


38473.3 


38027.8 






200 


63896.8 


62443.8 


61720.6 


50 


40 


100 


14368.5 


14038.7 


13874.6 






150 


38174.4 


37298.1 


36862.1 






200 


61943.9 


60521. 6 


59814.1 
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TABLE 4 



Theoretical Finned Condenser Heat Transfer 
Rate Versus Number of Axial Increments 

= 5000 BTU/hr-ft^-deg F At 3000 RPM 



Fin Half 


No . of 




No 


1 . of Axial 


, Increments 


Angle 


Fins 


S 


25 


50 


100 


10 


276 


100 


45685.4 


44652.3 


44138.1 






150 


120807.7 


118075.1 


116715.1 






200 


195703.4 


191270.2 


189066.4 


30 


84 


100 


42327.3 


41377.0 


40904.1 






150 


111871. 5 


109356.8 


108106.0 






200 


181198.5 


177121.2 


175092.6 


50 


40 


100 


39287.1 


38402.4 


37962.4 






150 


103724.0 


101386.0 


100222.1 






200 


167950.6 


164155. 9 


162268.8 
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t = 0.03125 in. 

b = 0.025 in . 

R = 0.762 in. 
o 

(}) = 



Figure 10. Condenser Geometry Considered in 
Two-Dimensional Analysis with 
3 Finite Elements 
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two equal parts, as shown in Figure 11, the geometry 
becomes 4 elements and 23 nodes. The computer program 
was modified again to solve this case, using the same 
boundary conditions and assumptions as before. Results 
for this case are also tabulated in Tables 5 through 10. 

C. DISCUSSION OF NUMERICAL RESULTS 

Computer results for the finned copper condenser were 
generated for three outside convective coefficients, three 
fin angles and two rotation rates. The working fluid used 
throughout the analysis was water. Tables 5 through 10 
compare the results of using 2, 3 and 4 elements, with those 
of Schafer [15 ] . 

Using the Finite Element solution, the heat transfer 

rates change less than 10 percent, as the number of elements 

change, so the heat transfer rates in the finned condenser 

can be approximated by 2 elements, and the error in the 

calculation of the heat transfer rates at the fin apex can 

be neglected. In comparing the smooth condenser to the 

finned condenser, the most dramatic improvement occurs 

for .a = 50 degrees with 40 fins. It is clear in this 

situation that Q. ./q ,, is far greater than 

4 Element smooth ’ 

^4 Element'^'^smooth * 

The heat transfer rate was found to vary slightly with 
rotational speed in a finned condenser. This result agrees 
with the results obtained earlier by Schafer [15] . Schafer 
suggested that the invariance in heat transfer rate with 
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Figure 11. Condenser Geometry Considered in Two- 

Dimensional Analysis with 4 Finite Elements 



64 



4 






« 



If 



) 




' 

I 

I 



I 

I 

4 




the rotation rate was due to the cancelling effects of 

decreasing film thickness and decreasing temperature 

difference across the film thickness as rotation rate is 

increased. The influence of external heat transfer 

coefficient, is shown in Tables 5 through 10. A 

large external heat transfer coefficient (h ^ = 5000 
^ out 

2 

BTU/hr-ft -deg F) is seen to produce the largest heat 
transfer rates. 
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III. EXPERIMENTAL PROGRAM 



A. EXPERIMENTAL EQUIPMENT 

1 . Description of Equipment 

The equipment used to support the experimental 
program was the same equipment used by Loynes [18] and 
Wagenseil [19], Figures 12 and 13 depict the experimental 
apparatus. The evaporator was a copper cylinder with 
0-ring seals at each end. Heat input to the evaporator 
was provided by a resistance heating element wrapped 
around the evaporator . The power supply for the heater 
element was changed from load bank resistors to a solid 
state phase amplifier power controller to reduce the space 
and inconvenience of the load banks. The model PA-1-2490 
phase/amp is a single phase silicon controlled rectifier 
power amplifier designed to control and regulate A-C power. 

A six-volt battery was used as a power supply for the 
controller circuit. 

The condenser was surrounded by a sealed insulated 
box. Gar lock seals at each end of the condenser kept the 
cooling water from leaking out. The cooling water system, 
as shown schematically in Figure 14, was four supply tubes, 
0.25 inch diameter with 0.031 inch holes drilled every 
0.20 inch. The supply tubes ran the length of the conden- 
ser, parallel to the axis. Filtered and softened tap water 
was used as the cooling fluid. The inlet and outlet cooling 
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Figure 12. Cross Section of the Rotating Heat Pipe 
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Figure 13. Schematic Diagram of Heat Pipe Apparatus 
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Cooling water 
inlet 




Section A-A 

Figure 14. Impingement Cooling Schematic 
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water temperatures were measured along with the cooling 
water flow rate. The heat pipe was rotated in two bearings 
by a three phase variable speed motor with a V-belt drive. 

A small drop-oiler was installed over the main 
bearing housing to lubricate the main bearing. This unit 
consisted of a small reservoir with a needle valve for flow 
regulation. A sight glass permitted observation of the oil 
flow from the reservoir to the main bearing. A flow rate 
of one drop every 20-30 seconds was found to be adequate. 

The end of the evaporator opposite the condenser contained 
a viewport consisting of two glass disks with an air space 
between them to prevent fogging. A strobe light shined 
through this viewport and, when synchronized to the rota- 
tional speed, allowed observation of the boiling action in 
the fluid annulus and the condensate in the condenser. 

2 . Instrumentation 

The condenser wall was grooved out externally in 
nine places in Figures 15 and 16. The grooves were spaced 
40 degrees apart and had a depth of 0.010 inch. Welded 
thermocouple beads were soldered into the grooves with a 
standard lead-tin low-temperature solder. The leads were 
secured to the condenser by banding with twisted wire and 
were passed through the small condenser flange to the mercury 
slip-ring unit on the test stand. Two additional thermo- 
couples sheathed in 0.0625 inch stainless steel tubing 
penetrated the main flange. The sheathes protruded 2.0 inches 
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Figure 15. Thermocouple Positions For 1.46 Inch Diameter and 1.00 Inch 
Diameter Cylinder 
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Figure 16. Thermocouple Positions For 0.5 Degree Cone 
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into the evaporator to measure the vapor temperature and 
were soldered closed at the tip to maintain the system 
vacuum. The inlet and outlet temperatures of the cooling 
water were monitored by Kapton insulated copper-constantan 
thermocouples. All thermocouples were read directly on a 
Newport Model 267 digital pyrometer. Calibration was 
performed in a Rosemount constant temperature bath using a 
platinum resistance thermometer as a standard. A graph of 
thermocouple error as a function of thermocouple reading 
was plotted to be used to correct each experimental reading. 
Cooling water flow rate was measured with a standard rotameter. 
Accurate regulation of the flow rate was accomplished by 
adjusting the pressure regulator installed in the cooling 
line . 

B . EXPERIMENTAL PROCEDURES 

1 . Preparation of Condenser Wall 

Tucker [20] established the procedure for prepara- 
tion of the copper condenser for obtaining film condensation 
with water. This procedure was followed with the modifica- 
tions listed below: 

a) The heat pipe was tilted so that the evaporator 
was depressed slightly. 

b) The standard cleaning steps were followed by 
scrubbing with a bristle brush. By depressing 
the evaporator end, cleaning fluids were prevented 
from running back into the heater assembly. 
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c) Complete wetting of the condenser wall was 
observed after cleaning. 

2 . Filling and Venting Procedure 

Instead of using the filling procedure of previous 
studies, a simpler procedure was used. The condenser was 
depressed, to nearly the vertical, with the evaporator end 
raised up. The heat pipe was then filled with about 350 ml 
of working fluid. The 0-ring and glass viewports were 
installed. The bolts of the retainer ring were torqued 
to thirty inch-pounds in ten inch-pound steps to avoid 
cracking the glass window. This filling procedure was 
followed by a venting procedure, used by Tucker [20] , with 
the modifications listed below: 

a) Tilt the heat pipe with the evaporator end 
down, until the working fluid covers up the 
viewport . 

b) Rotate the heat pipe by hand so that the vent 
screw is on top. 

c) Slowly apply the power to the heater to about 
0.25 KW until the saturation temperature in 
the. heat pipe is at least 220 degrees Fahrenheit 
for the distilled water. 

d) Back out the vent screw allowing some working 
fluid vapor to escape. Control the flow of the 
vapor with the vent screw, being careful to keep 
the pressure in the heat pipe above atmospheric 
pressure . 
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e) Vent for a minimum of ten minutes for distilled 
water, to allow escaping steam to drive off 
air trapped within the system, leaving about 
250 ml of working fluid left over. 

f) Turn off the electrical power and reseat the 
vent screw. 

g) Return the heat pipe to the horizontal position, 
and turn on the cooling water to the condenser. 

3 . Run Procedure 

With filling and venting completed, the system was 
ready to be run. The standard run procedure was as follows: 

a) Open cooling water supply to the condenser and 
bearing cooling water lines. 

b) Start the cooling water supply pump, and adjust 
the rotameter to the desired flow rate. In this 
experiment, only the 50 percent flow rate was used. 

c) Oil the evaporator bearing. 

d) Energize and start the drive motor. Increase the 
speed to 1100 RPM to insure an annulus has formed 
in the evaporator. 

e) Set rotational speed to the desired value. 

f) Let the system reach steady state; take zero 
power thermocouple, cooling water flow rate and 
tachometer readings. 

g) Set power to heater to desired level and allow 
the system to come to steady state. Take 
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tachometer, cooling water flow rate, thermo- 
couple, voltage and ampere readings. 

h) Repeat step g) for each power level desired. 

i) Repeat step f) . 

j) Repeat steps e) through i) for each rotational 
speed desired. 

4 . Data Reduction 

Data were reduced by correcting the recorded thermo- 
couple readings to obtain a true reading in degrees Celsius. 
Correction was accomplished using thermocouple calibration 
curves. The Celsius readings were then converted to degrees 
Fahrenheit. Percentage flow rate was converted to pounds 
per hour using the rotameter calibration curve. Heat trans- 
fer for each run was calculated as the product of the cooling 
water mass flow rate, the difference in temperature between 
cooling water outlet and inlet and the specific heat of 
the cooling water using the basic equation: 

Q = m C .AT (40) 

P 

where 

m = cooling water mass flow rate, Ibm/hr 

C = specific heat of the cooling water, taken 
^ as one BTU/lbm-deg F 

AT = the temperature difference between the 
cooling water outlet and inlet, deg F. 
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To correct for heat generated by friction in the 
Garlock seals and by viscous dissipation of the cooling 
water, the average zero power heat transfer rate was sub- 
tracted from the measured heat transfer rate in order to 
give a corrected heat transfer rate. 
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IV. DISCUSSION OF EXPERIMENTAL RESULTS 



A. EXPERIMENTAL RESULTS 

The various experimental runs are shown in Table 11. 

All the heat pipes in this thesis were run at 700, 1400 and 
2800 RPM, except the 1.00 inch diameter cylindrical condenser 
which was run at 1800 and 2800 RPM for the fourth run and 
1400 and 2800 RPM for the fifth run due to erratic readings 
of most of the thermocouples at the low rotation speed. 

It was concluded that the problem existed in the slip-ring 
unit and was due to dirty mercury in the slip-ring. So 
after the fifth run, the slip-ring unit was cleaned by 
taking the mercury and using methanol to clean the inside. 
Nitrogen was blown into the unit to dry it out. New mercury 
was added prior to the sixth run. The plots of heat transfer 
rate versus saturation temperature are referred to as the 
performance curves. 

Figures 17, 18, and 19 show the results of different 
condensers. In all cases, the heat transfer rate increased 
with increasing RPM. In the cases of the rotating cylindri- 
cal condenser, the condensate flow is induced not by centri- 
fugal force but by a hydrostatic pressure gradient established 
in the condensate as a result of its variable film thickness. 
By increasing the diameter of the cylinder from 1.00 inch 
to 1.46 inches, the heat transfer capacity increases. For 
the 0.5 degree truncated cone, the heat transfer is higher 
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TABLE 11 



Summary of Experimental Runs 



Run 


Geometry 


RPM 


Test Fluid 


No. 








1 


Cylindrical 1.46 


700 


Water 




inch diameter 


1400 

2800 




2 


Cylindrical 1.46 


700 


Freon 113 




inch diameter 


1400 

2800 




3 


Cylindrical 1.46 


700 


Ethanol 




inch diameter 


1400 

2800 




4 


Cylindrical 1.00 


1800 


Water 




inch diameter 


2800 




5 


Cylindrical 1.00 


1400 


Water 




inch diameter 


2800 




6 


0 . 5 degree trun- 


700 


Water 




cated cone 


1400 

2800 




7 


0 . 5 degree trun- 


700 


Water 




cated cone 


1400 

2800 
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Figure 17. Heat Transfer Rate Versus Saturation Temperature 
for Cylindrical 1.00 Inch Diameter Condenser 
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Figure 18. Heat Transfer Rate Versus Saturation Temperature 
for Cylindrical 1.46 Inch Diameter Condenser 
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Figure 19. Heat Transfer Rate Versus Saturation Temperature 
for Truncated Cone Condenser with 0.5 Degree 
Half Angle 
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than the cylindrical condenser. This improvement is due 
to the centrifugal force which increases the axial force 
component on the condensate and thereby increases the con- 
densate flow rate. This increase in working fluid circula- 
tion simply means that more vapor can be condensed per unit 
time. The increase in centrifugal force also flattens the 
condensate film and thereby reduces the overall thermal 
resistance in the condenser. 

The wall temperature profiles as shown in Figures 20, 

21 and 22 were taken from the 1.46 inch diameter cylin- 
drical condenser at different RPM. Figure 23 was taken 
from the 0.5 degree truncated cone at 700 RPM. Distance 
zero corresponds to the condenser end nearest the evaporator. 
The temperatures would be expected to be nearly constant 
if the heat pipe were truly isothermal. For the cylindri- 
cal condenser, the highest temperature is at the evaporator 
end, and the lowest temperature is at about 0.25 inch from 
the evaporator end. For the truncated cone, the tempera- 
ture profile decreases almost linearly, and the highest 
temperature is at the evaporator end. This increase in 
temperature at the evaporator is probably due to conduction 
effects from the heater and the main bearing. 

B. ANALYSIS OF THE RESULTS 

To compare the theoretical and experimental results, 
a log-log plot of the mean Nusselt number versus the Sherwood 
number was made. To make this comparison of theoretical 
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Figure 20. Condenser Outside Surface Temperature Profile at 700 RPM for 
1.46 Inch Diameter Cylindrical Condenser 
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Figure 21. Condenser Outside Surface Temperature at 1400 RPM for 
1.46 Inch Diameter Cylindrical Condenser 
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Distance Along Condenser, inch 

Figure 22. Condenser Outside Surface Temperature Profile at 2800 RPM 
for 1,46 Inch Diameter Cylindrical Condenser 
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Distance along Condenser, inch 

Figure 23. Condenser Outside Surface Temperature Profile at 700 RPM 
for 0 . 5 Degree Truncated Cone 
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versus experimental results, the following simplifications 
were made ; 

1) The thermal conductivity of the wall was constant 
along the condenser, 

2) the inside wall temperature was found by using a 
direct application of Fourier's law, 

3) all fluid properties were evaluated at the 
saturation temperature, 

4) heat flux through the condenser was considered 
uniform, 

5) for the truncated cone, the surface area used the 
average radius between the entrance and the exit 
radius , and 

6) for h = Q/AAT, AT was taken as the difference 

m 

between the saturation temperature and the condenser 
inner wall temperature. 

The mean Nusselt number for a given heat transfer rate 
is given by: 



Nu = h L/k 
m m 



(Q/A AT)L 
k 



(41) 



where 



Nu = mean Nusselt number 
m 

2 

h = mean heat transfer coefficient, BTU/hr-ft -deg F 
m 



= the average surface area, ft' 
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k = thermal conductivity of the fluid, 
BTU/hr-ft-deg F 

L = condenser length, ft 

The mean Nusselt number was determined from the data, 

and this value was plotted against the corresponding Sherwood 

number. Figure 24 shows Nu versus Sherwood number for 

m 

the cylindrical condensers. For water with the 1.46 inch 
diameter cylindrical condenser, the experimental data 
agree with equation (11). The data using the 1.00 inch 
diameter cylindrical condenser, however, show that the Nu^ 
increases with decreasing Sh. This trend is the opposite 
to what was expected. The reason for this trend may be 
due to: (1) ripples in the condensate film due to higher 

vapor velocities or (2) a high uncertainty in the low 
power runs. For Freon 113 and ehtanol, with the 1.46 inch 
diameter cylindrical condenser, the experimental results 
are higher than but are parallel to equation (11) . The 
reason for this discrepancy is not clear at this time. It 
may be due to differences in vapor densities between water. 
Freon 113 and ethanol, leading to differences in non- 
condensible gas migration, or to differences in vapor 
velocity effects within the theoretical model. 

Figure 25 shows the mean Nusselt number versus the 
Sherwood number of the truncated cone to find out the best 
correlation between the theoretical and experimental results. 
The results of Dhir and Lienhard [8], and Roetzel [14] were 
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Figure 24. Nusselt Number Versus Sherwood Number for Cylindrical Condensers 
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Figure 25. Nusselt Number Versus Sherwood Number for Truncated Cone Condensers 
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chosen to be compared with the experimental results. For 
the one degree truncated cone, using Tucker's data [20] 
and Wagenseil's data [19], the experimental results for 
water fall between both analyses, whereas for ethanol and 
Freon 113, the experimental results are higher than both 
analyses. This may be due again to differences in vapor 
densities between water. Freon 113 and ethanol, leading to 
differences in non-condensible gas migration, or to differ- 
ences in vapor velocity effects within the theoretical 
model. For the 0.5 degree truncated cone with water as the 
working fluid, the experimental results fall lower than both 
analyses, but fall higher than the analysis of Nimmo and 
Leppert [7] for a straight cylinder. 

Figure 26 shows Nu^/Nu^ plotted versus the parameters 
of Roetzel's analysis. By using all experimental results 
taken to data at the Naval Postgraduate School, it appears 
that good agreement exists between Roetzel's analysis and 
experimental results. The differences of experimental 
results and Roetzel's for the cylindrical condensers are 
less than 15 percent. For the truncated cone condensers, 
the differences are less than 20 percent for water. For the 
truncated cone with ethanol and Freon 113, experimental 
results are higher than Roetzel's analysis by about 35 
percent. However, the agreement is reasonable and the 
analysis of Roetzel can be used to predict rotating heat 
pipe performance to within ± 20 percent for water. For Freon 
113 and ethanol the data lies above the theory by +35 percent. 
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Figure 26 . Comparison of Experimental Data with Theory of Roetzel [ 14 ] 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

For an internally finned condenser, the analysis of 
the condenser heat transfer problem using a two-dimensional 
wall conduction model and a parabolic temperature distribution 
along the fin can be approximated by the two-element Finite 
Element solution with twenty-five axial increments. 

For a smooth condenser : 

1. The theoretical work of Leppert and Nimmo [6,7] can 
be used to predict satisfactorily cylindrical heat 
pipe performance. 

2. The theoretical equation which best predicts overall 
heat pipe. performance is due to Roetzel [14] ; 

XT M n ^ 4. A 23/16 ^.23/80, 4/23 

Nu = Nu (1 + 1.924 tan Sh ] 

m 

00 

The above equation agrees with the experimental data 
of water to within ± 20 percent. The experimental 
data of Freon 113 and ethanol lie approximately 35 
percent above the theory. 

B . REC OMMENDAT IONS 

Continue to study the internally finned, rotating heat 
pipe by: 
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1. Re-examining the Finite Element analysis using the 
mirror geometry presented in Figure 7. 

2. Analyzing film condensation on a variety of fin 
shapes, and applying the solution of the two- 
dimensional model to determine the shape permitting 
optimum heat transfer rate. 

3. Manufacturing an internally finned condenser to 
obtain experimental data to compare with the 
theoretical analysis. 

Continue to study the smooth, rotating heat pipe by: 

1. Using a condenser section with a half cone angle 
greater than 1 degree to obtain additional experi- 
mental data to compare to Roetzel's analysis. 

2. Examining reasons why the Freon 113 and ethanol 
data are falling above the theory. 
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APPENDIX A 



CALIBRATION OF MEASUREMENT DEVICES 

1 . Calibration of the Thermocouples 

The thermocouples for both the cylindrical and the 

truncated cone were calibrated, after they had been soldered 

to the outside condenser wall, in a ROSEMOUNT variable 

temperature oil bath utilizing a platinum resistance 

thermometer as the standard of measurement. They were 

calibrated in five degree increments, from 10 to 90 degrees 

Celsius. Each thermocouple was wired to a channel on the 

Newport Model 267 Digital Pyrometer giving individual 

readings in degrees Celsius. For each thermocouple the 

difference T . - T. was calculated where; 

act 1 

T . = reading in degrees Celsius of the standard, 

3 . j 

and 

til 

T. = reading in degrees Celsius for the i 
^ thermocouple at the same temperature. 

The mean of these differences was determined and 
used as the data point for plotting, - “^j-gcorded' 

versus thermocouple reading in degrees Celsius at known 
temperature. The maximum difference between the mean 
reading and each individual reading was one degree Celsius. 

2 . Calibration of the Rotameter 

The cooling water flow rotameter was calibrated 
using a Toledo scale, a large container and an electric 
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timer. Calibration proceeded in 5 percent flow increments 
from 15 to 90 percent flow rate. Water was collected for 
as long as 300 seconds at the lower flow rates and for a 
minimum of 140 seconds at the maximum flow rate. Each flow 
rate was measured two times. The flow rate (Ibm/hr) was 
plotted versus percentage flow through the rotameter. This 
curve was used to determine mass flow rate. 
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APPENDIX B 



UNCERTAINTY ANALYSIS 

The method of Kline and McClintock [21] was used to 
estimate the uncertainties in the experimental heat transfer 
rate. The heat transfer of the heat pipe is given by: 



^total 


= Q. = m AT, 

t p t 


m 

^frictional 


, = = m C AT_c 

loss f P f 


^corrected 


II 

o 

ft 

1 

o 



where 



II 

a 


total heat transferred to cooling water, 
BTU/hr 


lO 

II 


heat transferred to cooling water due to 
friction, BTU/hr 


Qc = 


corrected heat transfer rate, BTU/hr 


• 

m = 


mass flow rate of cooling water, Ibm/hr 


C 

P 


specific heat of cooling water, BTU/lbm-deg F 


t> 

ft 

II 


difference of cooling water outlet and 
inlet temperatures , deg F 


II 

<3 


difference of cooling water outlet and 
inlet temperatures due to heat of friction, 
deg F 
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The uncertainties of these quantities are designated: 




W , and 

m' C At 

P 



The fractional uncertainties are given by 





^"^ATf> 2 1/2 
AT^ ’ ' 



(B-1) 



^t m ^^t 



(B-2) 



The uncertainty for the corrected heat transfer rate is 
given by : 






(B-3) 



The values and uncertainties of the measured quantities 
were taken to be: 



Cp = 1.000 BTU/lbm-deg F 



= ±0.002 BTU/lbm-deg F 

P 



n 



= 50 percent of maximum flow rate = 1790 Ibm/hr 



W* = ± 1 percent 

m 



w w 

T. ' T , = ±0.18 deg F 

in out ^ 



W = (W^ + 

AT ^ T. -^T 

in out 
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Using the 1.460 inch diameter cylindrical condenser run 
at 700 RPM, the uncertainty limits are: 



T (deg f) Q (BTU/hr) UNCERTAINTY, W 

sc 



110.6 


1611.0 


± 634.6 


± 39.4% 


133.8 


3222.0 


± 637.7 


± 19.3% 


138.2 


3222.0 


± 641.4 


± 19.9% 


168.4 


5155.2 


± 643.6 


± 12.5% 


198.1 


7088.4 


± 651.8 


± 9.2% 


218.9 


7732.8 


± 655.0 


± 3.5% 


222.8 


3377.2 


± 658.4 


± 7.8% 
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